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Summary: 

In this paper we present a martingale related to the exit measures of super Brow- 
nian motion. By changing measure with this martingale in the canonical way we 
have a new process associated with the conditioned exit measure. This measure 
is shown to be identical to a measure generated by a non-homogeneous branching 
particle system with immigration of mass. An application is given to the prob- 
lem of conditioning the exit measure to hit a number of specified points on the 
boundary of a domain. The results are similar in flavor to the "immortal particle" 
picture of conditioned super Brownian motion but more general, as the change 
of measure is given by a martingale which need not arise from a single harmonic 
function. 
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1. Introduction 

One may think about super Brownian motion heuristically as a measure-valued process 
induced by a branching particle system in which each particle diffuses as a Brownian motion 
and critical branching is performed. With this image in mind, Dynkin's exit measure for 
super Brownian motion from a domain D is a measure on the boundary of D } supported on 
the set of points where the particles first exit the domain. Le Gall proved in [19] that the 
exit measure will "hit" a fixed point on the boundary of D } with positive probability, only 
in dimensions 1 and 2. In higher dimensions, a fixed point on the boundary will be hit by 
an exiting particle with probability 0. One may ask, if we condition the process to hit that 
point, what do the paths of the particles that hit the point look like? A naive guess is to say 
they are Brownian particles which are conditioned to hit the point, and so should look like 
/^-transforms of Brownian motion. In dimension 2, where no conditioning is necessary, this 
is essentially correct. Although the first hit, as chronicled by the Brownian snake, is not an 
/^-process, the typical path is. This has been exploited by Abraham [1]. 
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A rough heuristic shows that one's intuition should be true in dimension d > 3. Let z be 
a fixed point on dD } and for each e > define A e = A e (z) = dD D B(z } e) to be a ball on 
the boundary of D centered at z with radius e. Consider the "probability" u e (x) that the 
exit measure will charge the ball A e , under the excursion measure for the Brownian snake 
started from x. It is known (cf. Abraham and LeGall [2]) that if D is regular enough, then 
u f will solve the equation Au = 4u 2 , with boundary condition that is infinite on A e and on 
its complement. Fix some point x £ D. Then the ratio of functions v f (x) = u e (x) / u e (x ) 
will solve the equation Av e = 4v^u e (x ). As e — > 0, we will have that u f — > 0, and thus v f 
should converge to a harmonic function. The boundary conditions suggest that the limit 
function is the harmonic function associated with Brownian motion conditioned to exit D 
at z - that is the Martin kernel K Xo (-,z). 

When this heuristic is formally inserted into the known formula for the Laplace transform, 
one finds that the exit measure conditioned to hit the point z will look like an /i-path 
transform of the exit measure. Let D be a bounded domain in IR d , and let Dk C D be 
an increasing sequence of smooth subdomains with D = UkDk- Let X Dk denote the exit 
measure from under an excursion measure N X} let (•, •} denote integration and let ^ be a 
test function defined on dD^. Then formally, 

limN* (exp(-(X D *,<£)) | X D (A t ) > 0) = N Xo ((e X p(-(X D \cf>))(X D \ K Xo (-, z))). 

t— >o 

In this case the "harmonic" function is now H(ji) = (/i, K Xo (- } z)). As is usual with /i-path 
transforms, the change of measure is non-singular only for events prior to the exit time. 
Hence the appearance of X Dk in the statement. 

In addition, there is a concrete representation of this /^-transformed process. It is given 
in terms of an "immortal particle", or "backbone", consisting of a Brownian particle condi- 
tioned to exit D at z. The conditioned super Brownian motion is obtained by allowing this 
particle to throw off mass in a uniform Poisson manner. This mass then evolves as with an 
(unconditioned) super Brownian motion. 

Such /^-transforms have been studied by Roelly-Coppoletta and Roualt [24], Evans and 
Perkins [16], and Overbeck [22], while the immortal particle representation originates with 
Evans [14]. For example, in [16] and in [24], super Brownian motion conditioned on non- 
extinction is shown to be described in terms of an /i-path transform, where the harmonic 
function is now H(]jl) = (/i, 1). That is, if X t denotes the super Brownian motion in IR d at 
time /, then conditioning on non-extinction produces a Laplace functional of the form 

%(eM-{XtA)) = n x ((e W {-{XtA)){Xul)). 

The immortal particle representation in this case was obtained in [14]. For h a harmonic 
function, [22] defines the transform by H(ji) = The immortal particle representation 

in that case was obtained in Overbeck [23]. 

Our interest is in more general conditionings, in which one conditions, not just upon one, 
but upon several specified boundary points being hit. In particular, we will show that if one 
fixes n distinct points on the boundary of a bounded smooth domain D } and one conditions 
the exit measure to hit all of them, then the conditioned exit measure can be described in 
terms of an explicit martingale change of measure, or /^-transform, the transforms of the 
preceding paragraph being a special case. We will then give a backbone representation of 
this /^-transform, where now the backbone is a tree with n leaves, each terminating at one of 
the n boundary points. Our main purpose is to explicitly describe the evolution of this tree. 
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We note that in Etheridge [13], a different conditioning gives rise to a particle repre- 
sentation, where the backbone is also a tree. In Evans and O'Connell [15], a supercritical 
superprocess is given a particle representation in terms of a tree backbone throwing off sub- 
critically evolving mass. In Dembo and Zeitouni [9], a large deviation result with some 
similarities to our problem gives rise to yet another tree structure. 

We also note that steps have been taken towards a version of the Martin boundary, in 
the context of conditionings involving a single target (see Dynkin and Kuznetsov [12]). We 
conjecture that a more general Martin boundary theory would incorporate a modification of 
our basic conditionings, given in Section 6. 

To give a more precise formulation of the results, we first describe a general class of 
/^-transforms, which will include the ones we require as a particular case. Section 3 of the 
paper gives an analysis of this general class of transforms. To describe this class, we introduce 
some more notation. Let g > be a solution to the non-linear equation |Ait = 2u 2 in D. 
Let L± g be the generator of Brownian motion killed at rate 4g, and denote the law of the 
corresponding process by E 4a . Suppose we have n positive solutions in D to the linear 
equation L 4g v = 0, labeled v 1 , . . . } v n . Let U 49 be the potential operator for the generator 
L 4g in D and recursively define a family of functions v A indexed by nonempty subsets A of 
N = {1, . . . n}, as follows 

2E B a U 43 (v B v A \ B ) |A|>2.' 

For ^ A C N define 

M A = ]T ex V -{X\g)\[{X\v c ), 

aeV(A) ce<j 

where V(A) is the set of partitions of A. 

It turns out that {M A }k forms a family of martingales for the exit measure with cr-fields 
Tk generated by the superprocess paths before they exit D^. 

Theorem (Theorem 3.1). For each nonempty subset of A C N, M A is an T^-martingale 
under n x with N X (M A ) = v A (x). 

Changing measure by M A thus consistently defines the /^-transform 

M x (ex ? -(X k , ( f>)) = ^-^(expi-iX^^Mt). 

Such a transform appears when one conditions the exit measure to hit n distinct points 
on the boundary of the domain, as in the following result. 

Theorem (Theorem 5.6). Let D be a bounded Lipschitz domain in IR d ; d > 4, and let 
{ z i}?=i be n distinct points on the boundary of D. Set g = and let v l = K x (-,Zi) be the 
Martin kernel with base point, x and pole at Z{. Then 

■a 

llmN x (exp(-(X k ^)) | T[X D (A t (z i j) > o) = -^ N;c (exp(-(X fc , <j>))M» ). 

8 = 1 V ' 

Section 5 is largely devoted to proving this theorem. The arguments are based on an 
analysis of the asymptotics of small solutions to certain non-linear PDEs. Many of the basic 
estimates draw on the work of Abraham and LeGall [2]. In Section 6, we give a related 
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conjecture concerning conditioning super Brownian motion to hit the boundary in precisely 
n points. 

Having given an analytic formula for the conditioned exit measures, in terms of the mar- 
tingale change of measure by Mj^ , we turn in Section 4 to generalizing the immortal particle 
picture and representing dynamically in terms of a backbone shedding mass. This rep- 
resentation works for the general class of processes described above, with the effect of g 
being to prune off mass that would otherwise reach dD. To summarize, in Theorem 4.2 we 
construct a backbone, which consists of a branching process that eventually branches into 
n particles. In the example above these individual particles move as suitable transforms of 
Brownian motion. The particles in some sense keep track of a number of different exit points. 
Eventually there are n particles, each of which moves as a Brownian motion conditioned to 
exit D at a certain point. On top of this backbone, mass is immigrated in a Poisson manner 
and evolves as an unconditioned (but possibly pruned) super Brownian motion, to give us 
exit measures on {dD^}. In the case there is just one particle to hit, then the backbone is 
simply an /^-transform of Brownian motion. The arguments are based on an extended Palm 
formula for super Brownian motion. 

Conditioning the exit measure to hit a point z only produces a transform of the above 
type in dimensions high enough that the probability of such an event is actually zero. In 
dimension 2, any boundary point is hit with positive probability, and one must seek for 
other types of backbone representations for the conditioned exit measures. We give such a 
representation (among others) in an accompanying paper, [26]. 

For clarity, we will carry out the proof of the backbone representation in a simplified 
setting, avoiding the use of historical superprocesses. In fact, the arguments given would 
yield the stronger version of the representation, at the cost of a more baroque notation. In 
section 7 we will sketch the modifications needed to derive the strengthened results. 

Section 2 gathers together some basic facts which are used in the rest of the paper. 



2.1. Notation. For a set A, let \A\ denote its cardinality, and let V(A) denote the collection 
of partitions of A. Choose some arbitrary linear order -< on the set of finite subsets of the 
integers. For A such a finite subset, and a £ V(A), let cr(j') be the jth element of a in this 
order. Thus for example, 



We will switch between these notations according to which seems clearer. 

2.2. Set facts. We begin with some lemmas about sets. We will use the convention that a 
sum over an empty set is 0. 

Lemma 2.1. Let A C B C C be subsets of {1, 2, . . . n}. Then 



2. Preliminaries 



(7 



Y[(X k ,v c ) = \[{X\v a{ ^). 



C'Ea j = l 





A=C- 



ACBCC 



Proof. 



i)W = (_i)i*i £ (-i)W 



ACBCC 



BCC\A 
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f(-l)l c l C = A 

\{-l)\ A \{l-l)\ c \ A \ C\A^0. 



□ 



The following is immediate. 



Lemma 2.2. Let A be finite, and let Wi G R for i G A. Then 



n(i--o=i+E(- i ) |ci (n^) 



ieA cca iec 

He 



In this paper we use the letter K to denote a generic non-trivial constant whose particular 
value may vary from line to line. If it is important, explicit dependencies on other values 
will be specified. 

2.3. Facts about conditioned diffusions. First we recall some formulae for conditioned 
Brownian motion. 

Let B be (/-dimensional Brownian motion started from x, under a probability measure P x . 
Write To for the first exit time of B from D. 

Let g : D — > [0, oo) be bounded on compact subsets of D } and set 



Let ^ t be a process which, under a probability law PJ?, has the law of a diffusion with 
generator L g started at x and killed upon leaving D. In other words, £ is a Brownian motion 
on D } killed at rate g. Write ( for the lifetime of £. Then 



Let U 9 f(x) = Jq 00 P 9 ( / {£>t})c/t be the potential operator for L g . If g = we write U 
for U 9 . If < u is L 3 -superharmonic, then the law of the u-transform of £ is determined by 
the formula 



for $(£) G c{^ s ; 5 < t}. Assuming that < u < oo on D } this defines a diffusion on D. If u 
is L 3 -harmonic, then it dies only upon reaching dD. If for / > 0, u = U 3 f (that is, u is a 
potential) then it dies in the interior of D } and in fact P£' u satisfies 



L g = 2 A - 9- 








where is the process £ killed at time t. 
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2.4. Facts about the Brownian snake. Next we recall some useful facts about the Brow- 
nian snake. Refer to [7] or [10] for a general introduction to superprocesses. 

The Brownian snake is a path- valued process, devised by Le Gall as a means to construct 
super Brownian motion without limiting procedures. Refer to [18] or [20] for the construction. 

We use the standard notation (W S} ( S ) for the Brownian snake, and for the excursion 
measure of the Brownian snake starting from the trivial path (w } £), ( = 0, w(0) = x. Note 
that W s (-) is constant on [£ s , oo), and (. has the distribution of a Brownian excursion under 

Super Brownian motion X t is defined as 

(X t ,<f>) = j <f>(W s (t))dL t (s), 

where L t is the local time of (. at level t. Dynkin [11] introduced the exit measure X D 
associated with X t . We follow Le Gall's snake-based definition of X D (see [20]) as 

D i\ _ f iijjT (/■ \ \ jtD 



(X»,<j>) = J ^Ws{Cs))dL»(s), 

where L D ( ) is an appropriate local time for W s (( s ) on dD. (In section 7 we define these in 
a historical sense.) 

We denote the range of the Brownian snake by 1Z(W) = {W s (t) : < s < cr, < t < ( s } 
and the range inside D by Tl D {W) = {W s (t) : < s < cr,0 < t < t d (W s ) A ( s }. There is an 
obvious inclusion between the range inside D and the exit measures, given by 

{(X D , \ A ) > 0} C {R D (W) n A + 0}. 

We make use of the following facts about the Brownian snake, which are contained in 
Le Gall [20]. 

A useful first-moment calculation is: (cf. [20] Proposition 3.3) 

N x ((X D , ( f>)) = E x (cf>(B TD )). 

The following is an immediate consequence of Theorem 4.2 and its corollary in [20]. 

Lemma 2.3. Let g be a solution to Ag = 4g 2 in D , and let {-Dfc} be an increasing sequence 
of smooth subdomains of D. Then for each k, 

n x (l-exp-(X D *,g))=g(x). 

Remark 2.4. We make use of a generalization of this, replacing Brownian motion by Brow- 
nian motion killed at a certain rate. 

Let Th = ^D k be the cr-field of events determined by the superprocess killed upon exiting 
Dk- See [11] for a formal definition, or refer to section 7, where we give a definition in terms 
of the historical superprocess. 

Dynkin introduced a Markov property for the exit measures in [11]. In our context, the 
Markov property is established in [21]. The next result gives it in the form we will use it: 

Lemma 2.5. 

N*(exp ~{X D } <f>) | T k ) = exp -{X D « , N.(l - exp -{X D } </>))). 
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We use the following notation, where B s denotes a path in D whose definition will be clear 
from the context: 



e D 



B)=exp- I ds 4N Bt (l - ej). 
Jo 



The Palm formula for the Brownian snake takes the form: (cf. [20], Proposition 4.1) 

N,«X D ,</>) e J) = E x ^{B TD )K D {e°)). (2.3) 

We will make use of the following extension to the basic Palm formula. See [8] for a general 
discussion of this type of Palm formula. 

Lemma 2.6. Let N = {l,2,...n} ; n > 2. Let D be a domain, and let B be a Brownian 
motion in D with exit time r. Let {ipi} be a family of measurable functions. Then 



n x (e^l[{x D }1 p t }) 

ieN 

= \T, e x {a f dtM t (e^)n Bt (e^\[(x D ^))nB t (^ II ( xD ^)))- 

MCN ^° ieM ieN\M 

0,JV^M 

Proof. Observe first that, by monotone convergence, it suffices to prove the result in the case 
that D and the ipi are bounded, and (p is bounded away from 0. So assume this. 
Let N* = {2, 3 . . . , n}. Let D\ denote the derivative in A. Then 

M^l[(X D ,^)) = n x ((X D ^ 1 )e c , H(X D ,^)) (2.4) 

ieN ieN* 

= (-lT-'D^ ■ ■ ■ D Xn N x Vi>e(</> + X ^ U 2 =-=a„=o 

= (-l)"" 1 ^ • • • D Xn E x ^MB T )M r (e(<t> + X ^)fj U 2 =-=a„=o (2.5) 

= ]T E x U l (Br)M T (e 4> )\[U j dtN Bt (e^ ]J(X D ^ t )) 

Line (2.5) follows from the basic Palm formula. This last line uses the notation V(A) 
for all the partitions of a set A, and follows from a simple induction argument for the 
derivatives. Because D is bounded, differentiation under the integral sign is easy to justify 
using dominated convergence, once we establish that R^e^ JT^X" , ?/>;}) is bounded. But 
since the ipi are bounded, and cp is bounded from 0, in fact riiv*^ 13 ' V^'} 1S itself bounded, 
so that this follows from (2.3) and the boundedness of ip\. 

We use the notations abc ... to denote all the elements abc . . . but the 6th one. Let Q t be 
the filtration of B t . Then by the Markov property for Brownian motion, one has that 

n x {^l[(x D ,^)) 

ieN 
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E^MBrW^) 

aeV(N*) 



o Jo 



dh ■ ■ ■ dt wl Y[4N B{tj) (e^ Y[ (X D ^ t )) 



3 = 1 



ie<r(j) 



E x (MBr)KM52 f dt k An B[tk) (e^\[{X D ^ l )) 

a£V(N*) k=l J ° i£a(k) 

••• / dU---di k ---dt H [j4N B ( t .)(e^ \[{X D ^i)) 

Jtk 3 + k ie<r(j) 

m ( r 

J2 E -\ 4 dt k n B{tk) (e^ ]J{X D ^ t ))E x [ip 1 (B T )A r T (e^) 

a£V(N*) k=l \ J ° i£a(k) 

x j ■■ [ dh ■ ■ ■ dt k ■ ■ ■ dt W] l[4N B{tj) (e^ J] (X D ^ t )) | g tk U 

Jtk Jtk 3 + k ie<r(j) J 

E x Uj dtn Bt (e^ J](X D ,V 8 ))M(e^) 

MCN* \ JO icM 



0^MCJV 



x E Bt (MBr)M T (e^)l[(4 l' T dtn Bt (e^ J] (X D ',^») 



0^MCJV* 



iGJV\M 



MCN J° iC-M icN\M 

0,JV^M 

Note that in order to rewrite the set of partitions of N* } we designate one of the elements 
of a given partition as M, so that the rest forms a partition of N \ M . This isn't one-to- 
one, which is taken into account in the counting. The last factor of 1/2 occurs because the 
summation in the last line counts everything twice. □ 

We use the extended Palm formula to show an exponential bound on the moments of the 
exit measure. 

Lemma 2.7. Let D be a domain in W d satisfying sup xeD E x (td) < oo, where is the exit 
time from D for Brownian motion. Then there exists A > such that 

swpN x (exp\(X D ,l) - 1) < oo. 

xED 

Remark 2.8. A bounded domain D in IR d will satisfy sup D E x (td) < oo. 

Proof. It is enough to show that N x ((X D ,l) n ) < n\M n for some M < oo. Our proof 
follows part of the proof of Lemma 3.1 in Serlet [27]. Set c n = sup^ N x ((X D ,l) n ). Then 
C\ = sup D E x (td) < oo by assumption. By Lemma 2.6, with </> = we have the immediate 
recursion relation 

Ci = Ci 
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c n < 2ci ^2 (^J c n-3 C 3, n>2. 



Let a n+ i = Ci(2c\) n (2n)\ / n\ , a\ = c\. Then one can easily verify, by induction, that the 
sequence a n satisfies the combinatorial identity 



n-l 

n 



n — i / 

a n = 2c i y^{ 



Since C\ = a\ < oo, we get by induction that c n < a n . Stirling's formula applied to a n shows 
that some M < oo exists. □ 



3. A GENERAL CLASS OF MARTINGALES 

We now present a general class of martingales related to the exit measure. 

Let D be a domain in IR d , d > 1 and let g > be a solution in D to the non-linear equation 

^Au = 2u 2 . (3.1) 

In particular, g satisfies the conclusion of Lemma 2.3. Readers may feel free to take g = 0. 
In fact, the only application given in this paper requiring a non-zero g is that of Section 6. 
Another reason we opt to work with a general g is for consistency with [26], in which such 
g } s play an essential role. 

Recall that E 43 denotes the law of Brownian motion killed at rate 4g, and that L 4g is the 
generator of this process. That is 

L 4 gW = ~ ^9 W - 

Suppose we have n positive solutions in D to the linear equation L 4g v = 0, labeled v 1 , . . . , v n . 
Recall that U 43 is the potential operator for the generator L 4g in D } and recursively define 
a family of functions v A , for ^ A C N = {1, . . . n}, as follows: 



yA = ^ 2 Z bca U 43 (v B v A \ B ) \A\ > 2. ( 3 - 2 ) 

Note that v A is either finite everywhere on D, or v A = oo. 

Let Dk fT D be an increasing sequence of bounded, smooth subdomains and for each k let 
X k be the exit measure from Dk, and = exp — (X k , cj)). Let denote the exit time of a 
path from Dk- For |/AC]V define 

M£= eM-(X k ,g))l[(X k ,v c ). (3.3) 

aeV(A) ce<j 

It turns out that {M A }k forms a family of martingales for the exit measure. 

Theorem 3.1. For each nonempty subset of A C N, M A is an Tk-martingale under N x , 
provided v A < oo. 
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A = {i}, 



Remark 3.2. Because M A is a martingale we can define a martingale change of measure in 
a canonical way as follows. For a .^-measurable function, we set 

f 



v A (x) 



(3.4) 



Remark 3.3. Letting g = be the trivial solution to (3.1) and v = 1 (or v = h with h 
harmonic) we recover the well-known fact that (X k , 1), (or (X k ,h)) is a martingale. 

Proof. We first establish the following lemma 

Lemma 3.4. If v A < oo ; then N X (M A ) = v A (x). 

Proof. The proof relies on the extended Palm formula. We induct on the size of A. First 
consider the case \A\ = 1. 

N X (M A ) = N x (exp(-(X\g))(X\v A )) 

f' T k 

= E x (v A (C Tk )exp- 4% (1 - exp-(X\g))) (Palm formula) 

Jo 

f' T k 

= E x (v A (^ Tk )exp- ds4g(Q) (Lemma (2.3)) 

Jo 

= E^(v A (U),C>r k ) = v A (x). 
Next, suppose \A\ = n and that the lemma holds for all smaller sets. Then 
v A (x) = 2 U 49 (v B v A \ B )(x) 

BCA 
0<\B\<\A\ 

= E 4 j(v A (U)X< n)+2 ]T E x [ f dt exp(- f ds4g(Q)v B (C t )v A \ B (C t )). 

BCA ^° ^° 

0<\B\<\A\ 

Denote the first of the above terms by 7, and the remainder by 77. As in the case \A\ = 1, 
I = N x (e k (X k ,v A )). Further 



77 = 2 ]T E x (rdtAf t (e k g )Nt t (M»)Nt t (MP B )) 



BCA 
0<\B\<\A\ 



2 £ E x 

BCA 
0<\B\<\A\ 



(induction) 



dtM t (e k ) ]T ^e k g X[(X k 

veV(B) 



i=l 



(Definition of M A ) 



e %( e "n( xfc ^ 7(j) )) 

1CV{A\B) 3=1 



E E E ^. 

BCA (7e 'P(_B) 7e 'P( y l\_B) 
0<|B|<|A| 



N, 



M e 3 



i=l 



3 = 1 
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E E ^ 

aeV(A)CC{l,...\a\} 
0<|C|<k| 



dtM t {i 



iec 



]T N x (e k g Y[(X k ,v" 



(Palm formula) 



<reV(A) 

kl>i 



i=l 



Thus 



v A ix) = I + II= ^ N x (e k g l[{X k ,v°U)) = N x (M£). 



i=l 



□ 



Next we establish another lemma which helps us identify N X (M^ +1 \ ft). 
Lemma 3.5. 

i = e k g e ri<^> N -( e r n (x k+i ,v° u) ))) 

aeV{A) 1 eV{{l,...,\a\}) i=l jei(i) 

Proof. This lemma follows from the strong Markov property satisfied by the exit measures 
(Lemma 2.5) and differentiation. Differentiation under the integral sign can be justified as 
before. 



n x (M* +1 1 T k 



E (-1 PD M ■ ■ ■ D Xw] n x (exp -(X k+1 ,g + ]T A^«> | ft 

<tEV(A) i=l 

(-l) H D Xl ...D Xwi e X p-(X k , 

veV(A) 

M 



k J |Ai = -=A H =0 



(3.5) 



N. (l - exp - (X k+1 ,g + J2 X ^ {t 
exp-(X k ,K(l-exp-(X k+1 ,g))) 



|Ai = -=A| CT |=0 



<reV(A) 



(3.6) 



1 ev({i,...,\a\}) i=i je-y(i) 

= < E E n<**,N.( e ; +i n<^ +1 ^ (J) »- 

aeP(A)7eP({i,...,|a|}) i=i jei(i) 

Line (3.5) follows from the Markov property of the exit measures, and (3.6) is from the chain 
rule of calculus. □ 
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Finally, we prove the theorem. We have, by rearranging the following sum, that 



<reV(A) i=l 
a 

=< e n< x *> N -«S)> 

aeV(A) i=l 



(3.7) 



=< e n e (^,N.( e j +i n(^ +i ^ 7(j) ))) 

aeV(A) i = l j = l 

a 

=< e e n< x *> N -( e J +i n (x k+i ,v^))) (3.8) 

7G-P(A) aeP({l,- ,| 7 |}) i=l ie<r{i) 

= N(M A +1 I ^*). (3-9) 

Line (3.7) follows from Lemma 3.4, whereas line (3.8) follows from rearranging the terms 
in the sums. Line (3.9) is from Lemma 3.5. □ 

4. A BRANCHING PARTICLE DESCRIPTION 

We now show that changing measure via is equivalent to taking a branching particle 
backbone process with immigration of mass along the paths of the particles. To formulate 
this, we define two measures, and N^. Recall that the former is defined by 

dWL 1 

K- — 



so that 



M^exp 



v N (x 



-Mi 



N 



-N,(exp(-(A te ,</>»Mf) 



v N (x 



To define N X} we first construct the branching process which will form the backbone of 
our new exit measures. Let D, Dk, N, g and {v A } be as before. By (3.2) and the linearity of 
U 49 , we have that v A is an L 43 -potential for \A\ > 2 (provided it is finite). For \A\ = 1 it is 
iv4 3 -harmonic. The branching-particle process we desire will start with a single u^-particle, 
that is, a particle moving as a u^-transform of the process with generator L 4g . If |iV| = 1, 
then because v N is L^-harmonic, the particle dies on the boundary of D. If |iV| > 2 then v N 
is an iv4 3 -potential and the particle dies in the interior of D. At its death, it splits into two 
independent particles, a u^-particle and a v N \ A particle, with position-dependent probability 
given by 

(v A v N \ A )(y) 



p(A,N)(y) 



E 



(v B v N \ B )(y)' 



4.1 



The v A particle tracks the functions jV}; ey i while the 



N\A 



particle tracks the remaining 



functions. This pattern then repeats for each new particle. For example, the v A particle dies 
on the boundary of D if A is a singleton, and otherwise it dies in the interior of D } giving 
birth to two new, independent particles as before. Since the lifetime of each particle is finite 
almost surely, eventually all of the particles will have died out. In total there will be 2n — 1 
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particles of which n die on the boundary of D. The remaining n — 1 will die in the interior 
of D. 

Let n t denote the number of particles alive at time t. Label them with 1 < i < n t} and 
for each one set x 8 (s),0 < s < t to be the history (including the ancestors' history) of the 
individual particle up until time t. Define measure- valued branching processes as follows 

n t n t 

T t (dx) = 'Y^5 x i( t )(dx), T k (dx) = ^\ Tk{xJ)>t 8 XJ{t) (dx). (4.2) 

8 = 1 8=1 

The process T k puts a mass at each particle alive at time t which hasn't already exited D k . 
Without comment, these processes will be referred to in terms of the underlying particles 
although strictly speaking they are measures. Let Q x denote the law of T t . 

We now specify how the mass T throws off evolves. The following calculation shows that 
exp — (X k ,g) is an jF fc -martingale (although it has an infinite moment under N x ), and follows 
by the Markov property and Lemma 2.3. 

n x (exp-(X k ,g) | T k _ x ) = e W -{X k -\n.(l-e W -{X k ,g)) 

= exp-(X k -\g). 

Thus we can define a consistent measure on the T k measurable sets by 

N*($ fc ) = n x (<S> k exp-(X k ,g)). (4.3) 

The extra factor prunes off mass that would otherwise get to dD k . Though we will not 
need it in the rest of the argument, the following is a justification for this interpretation of 
N^. It is essentially a special case of Dawson's Girsanov formula (see section 10.1.2 of [7]). 

Lemma 4.1. Let be the excursion law for the superprocess in D based on the generator 
L 4g . Then N E (1 - exp -(X k , </>)) = N E (1 - exp -(X k , </>)) for every </>>0. 

Proof. Let ip be the solution to 

L 4g u = 2u 2 in D k 
u = (f> on dD k . 

It is then easily checked that u = ip -\- g is the solution to 

1 

L 

2 



— Au = 2u 2 in D k 



u = (f) -\- g on dDk, 

Thus by Lemma 2.3, 

n x (l - exp -(X k , </>)) = n x ((1 - exp -(X k , </>)) exp -(X k , gf) 

= n x (l- exp -(X k , </> + g))- n x (l - exp -(X k , g)) 

= (i>(x) + g(x)) - g(x) = i>( x ) 

= k(l-exp-(X k ,<f ) )). 

This suffices. □ 
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One thinks of mass being created continuously along the backbone, but only at countably 
many times will it actually survive, even instantaneously. At each such time, the mass 
created evolves like a superprocess with "law" some N y} and, if it survives long enough, 
produces a contribution to the exit measure. More properly, given the backbone T fc , we 
form a Poisson random measure N k (d/i) with intensity dt j 4T k (dy)'Ny(X k £ o?/i). We 
then realize the exit measure under as Y k = j ^N k (d^). A standard calculation shows 
that 

/>oo 

%(exp-{Y k }( j ) }) = Q x (exp- / dt 4{T k } N.(l - exp -{X k } <£)))). (4.4) 

Jo 

Theorem 4.2. If v N < oo, then ML, (exp -(X k , </>)) = N E (exp ~(Y k , </>)) for each </>>0. 
Remark 4.3. In section 7 we will sketch a proof that in fact, = on Tk- 
Proof. Recall that 

M E (exp-(X fc ,^)) = -±-^N x (e k M? ). 

For A C N, a £ V(A) define m£ = e k g \{ B€(r {X k ,v B ), so that A# = Z<,ev(A)™k- Fot 
B C A define 

. . . . , B m } B = U^Bj, each B 3 £ a 

o" |g = < J 

10 if no such decomposition exists. 

That is, if a can be restricted to be a partition of B } then define it as such, otherwise set 
cf\b to be the empty set. 

We use £ for the canonical process with lifetime (. 

We first mention a consequence of the Palm formula. Let a £ V(N). Then 

cc{i,...,W\} iec 

i<|c|<kl 

x%(4 + , II <**X W » 

t€{l,...,k|}\C 



1 



2 ^ 



^ E x (4 / ^M(eJ +fl )N 6 (eX U )%(eK'^))- ( 4 - 5 ) 



The branching processes T naturally partitions the \N\ points in the following manner. 
We eventually have n particles that exit the domain D. Group together those that have a 
common ancestor that exits Dk- We will write T k ~ a for the relationship that the resulting 
partition is a. 



Lemma 4.4. On the event {T k ~ a} one h 



as 



oo 

k 



Q x (exj>- / <ft4(T*,N.(l - exp </>))); T fc ~ a) = —-® x ( e * m Z). (4.6) 

Jo v \ x ) 
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Proof. This will be established by induction on the size of |<r|. First we look at the case 
|cr| = 1, in other words, where there is no branching and a = {N}. Since is formed by a 
u^-particle which exits D k before dying we have 

/>oo 

Q x (exp- / <ft4(T*,N.(l - exp -<X fc , </>))), T fc ~ a) 
Jo 

= E^ vN (ex ? - r^4%(l -exp-(X fc , </>», ( > r k ) 
Jo 

'•Tk 



J—Et a (v N (ti Tk )exv (- jf k dtm Ct (exp-(X k ,g)(l -exp-(X fc , </>»)), ( > r k ) 



-J—Ejv N (Z Tlc )exp(- [" dt4g(£ t ))exp( [ * dt 4% (1 - exp -(X\g))) 
v \ x ) v Jo Jo 

f' T k \ 

xexp(-y c/t4%(l - exp-(X\(f> + g)))J 

-^^^E x (v N (( Tk )ex V (- / dt4g(£ t ))exp( dt4g(£ t )) 
v \ x ) v Jo Jo 



expf 



C /t4%(l-exp-(X fc , <f> + g)))) 



1 -E x (v N (U)K k (ei +g )) 



v N (x)~ 

1 N x (e k e k g {X k } v N }) (4.7) 



v N (x) 

1 



-N x (e k ml 



Line (4.7) follows from the Palm formula with a = {N}. 

To establish (4.6) for all a £ V(N), we assume the induction hypothesis is true for all 
A C N and all partitions of A of size j or smaller, j ' > 1. Now consider cr £ V(N) of size 
J + l- 

Since |cr| > 2, we must have that T branches before it exits D k - The first branch time of 
T is the lifetime ( of the u^-particle. 

Q x (exp (- dt 4(Tf, N.(l - exp -(X\ <f>))fj ; T k ~ a) 

= Q,(exp (- Qf' + <ff4(T*,N.(l - exp -<X fc , ; T fc ~ a) 

C 

/I "R\T (-^ I Vk 1\ \ \ A AT\ff 



E^ N (l (<Tk exp(- [ C dtm Ct (l-exp-(X\ cf>))) ]T p(A,iV)(6 
V Jo J A ^ AT 



ACN 



Q Sc (exp {- dt 4(T k , N.(l - exp -<X fc , ; T fc ~ aU) 



(4.8) 



CO 



Q Sc (exp(- / <ft4(T?,N.(l -exp-(X fc , </>»> ) ; T fc ~ a^)) 
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E 4 J> vN (l (<Tk e X p(- [ C dtm Ct (l-e X p-(X k , <j>))) ]T P(A,N)& 

\ JO / A^AT 



ACN 
l,A^N 



(4.9) 



]T E 4 J> vN (l (<Tk e X p(- 4%(l-exp 

A r- AT \ J 



ACN 



dt4N 5t (l-exp-{X\(j)}) 



ACN 
*U^0 



^N £ ,(eXhN £ ,(^)(2 E (» fl ^ fl )&))) 



(4.10) 

2 Y, ^Tx) Ex{ l " ds eXp (" / ^4flf(6))exp( /" C /t4%(l-exp-(X fc ,(/))) 



0,A^JV 

<^U^0 



exp (-£dtWt t (l -exp-(X fc , + % (eJm^)Nc. (ejm^)) 



(4.11) 



2 £ ^^(r da ^( e 5 + J%( e X u )%( e K'^)) 

u \ x ) Jo 



ACN 
0,A^JV 

<^U^0 

= N s (eX). (4-12) 

Line (4.8) follows from the strong Markov property for the branching process T applied 
at the first branch time, and the description of the branching process having independent 
offspring. The next line, (4.9) then follows from the induction hypothesis. The formula (2.2) 
for a conditioned u-processes gives line (4.10) and (4.11) follows since under E 4a the process 
is Brownian motion killed at rate 4g. Finally, (4.12) comes from (4.5). □ 

By summing (4.6) over a £ V(N) we get 

k tSj /1 „„„ / vk 



Q x (exp- / <ft4(T?,N.(l -exp-(X te , </>)))) 
Jo 

^ ' aeV{N) 
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Using the definition of N from (4.4) we get that 

N*(exp -(Y k , </>)) = M^exp -(X k , </>)). 

□ 

5. Conditioning the support of the exit measure to hit n points on dD 

In this section we investigate the exit measure when it is conditioned to charge n small 
balls on the boundary of D. In the limit as the radius of the balls tends to 0, the conditioned 
process converges to one given by a martingale change of measure as in section 3. 

Let N = {1,2, . . . , n} and let {zi}i e ^ be a finite set of distinct points on the boundary 
of a bounded, Lipschitz domain D. Using the notation B(x, e) for a ball centered at x with 
radius e, define the sets 

Bi = B( Zi ,e), 
A\ = B{ n dD. 

For a set A C N set Bf = VJaB\ and Af = U^A*. Furthermore, for x £ D define the 
functions 



u e (x ) = n x (2^{X u } l A ( Zt ,e)} > 0) = N*(Hit at least one of A(zi, e) while exiting) 

vf(x) = n x ( Y[(X D , lA(z„e)} > 0) = n x (Kit all the A(zi, e) while exiting). 

ieA 

It follows from Theorem 7.1 of [2] that 

Uuf = 2(uff 



Al , i z £ dD \ Bf 
km uf(y) = < A 



A simple relationship between the two functions is summarized in this lemma which is 
proved by an inclusion exclusion argument. 

Lemma 5.1. We have the following relationships: 

vf = - (-1)1 V 

Q^BCA 

We now show that the family of functions vf satisfies the following recursive relations. 

Theorem 5.2. For D a bounded Lipschitz domain, and A and vf as above, the functions 
vf satisfy 

(a) vf £ C*(D). 

(b) \^vf = -2 £ (-l) |SnC W. 

9^B,CCA 
BUC=A 
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2((2«f + (-l)IV)^- E (-l) |SnC W 

|A|>2 




(C) y J£an V ?M = )r \A\ = 1. 

[ [oo z £ A A 

Proof. For convenience, drop the e subscript. When |A| = 1 this is Theorem 7.1 of [2]. Let 
denote surface measure on dD, which exists by the Lipschitz condition. The condition on 
D in [2] is that it satisfy 

P)(A a U A A ) 
liminf / a , > 0, for all a £ A^. 

r^O 0(A«) 

For this problem, this is trivially true as the limit is 1. 

For \A\ > 2, we have by Lemma 5.1 and the relation (1/2) Au A = 2(u A ) 2 , that 



0^£CA 

"2 E E (-1) |C| - C ) 2 (5-1) 

_ 2 £ (-1)1*1 ]T (-l^l+I^V 3 

0^£CA 9^C,DCB 

0^C*,DCA CuDCBCA 

"2 E (-l) H+|CnD| ^ D (-l) H lcuD^ (5.2) 



0^C*,DCA 



-2(2 e (-i) |c| ^-(-i) |yl| (^) 2 ) -2 e (-I)'--" 

$^CCA C'UD=A 



2(2(-«^-(-l)W(^) 2 )-2 E 



^ CnD \v c v D . 



CUD=A 
%,A+C,D 



Here line (5.1) follows by Lemma 5.1, and line (5.2) by Lemma 2.1. 

Finally, we remark that v A (y) < u 8 (y) for i £ A. Fix z £ dD. For e sufficiently small, 
because of the Lipschitz assumption on D } z is not contained in B\ for some i £ A. Thus as 
y — > z we, have < u e A (y) < u*(y) ^ by the |A| = 1 case. 

□ 

Theorem 5.3. Let D be a bounded Lipschitz domain in IR d for d > 4. Let {^i, . . . £ n } 6e 
distinct points on the boundary of D, A C N = {l,2,...n} ; and x } y £ D. Let £7 6e £/ie 
potential operator for Brownian motion killed upon hitting the boundary of D, K® the Martin 
kernel for D. Then, the following limit exists 

lim V ^ V) = m A (y). 
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Furthermore, the limit satisfies 

A , , \K°(y,zi) A = {i}, 

*^ y ) = 2Efla« fl )(y) 1^1 > 2- 



Proof. Fix 8 > such that B(z t} 8) B(zj,8) = if i ^ j £ A. Set Df = D \ U A B(z t} e) } 
and r e = = t^a. 

The proof of this theorem relies on the following lemma. 

Lemma 5.4. Assume the conditions of Theorem 5.3, and fix an x £ D. There exist con- 
stants K < oo and e > such that for all e < e , A C N and y £ Df^ t 



n 



true 



7n particular, the expression on the left is bounded in y on any compact subset of D, uniformly 



in e. 



Proof (of lemma). We prove the lemma by induction on the size of A. First the case A = {i}. 
Set z = Zi, and v = v\ This part simply extends the argument of the proof of Theorem 3.1 
of [2] to the Lipschitz case. 

Since Av e = 4(u e ) 2 in D and v f vanishes on dD \ A(z, e) the Feynman-Kac formula gives 



v t (y) = P y (v e (B(T 2e j) exp(- / 2v e (B s )d. 



y ' 

/o 



is 



< P y \v £ (B(t 2£ )), B(r 2e ) £ dB(z, 2e)J . (5.3) 

But, from the definition of v f one has the bound 

v e (w) = n w (X D (B(z,e))>0) 
< N w (K D (W)ndD + 0) 

<cd(w,dD)- 2 . (5.4) 

Line (5.4) follows from comparing the hitting probability with the exit probability from a 
ball (see Proposition 2.3 of [20]). 

Combining this with (5.3) and the definition of yields 

v e (y) < Ke- 2 P y [B(r 2t ) £ dB{z,2tj). 

We now show that for y £ Df f , 

P y {B T2e £ dB(z, 2e)) < KP y (B TD £ A(z, e)). (5.5) 

This leaves us with the bound (for d > 3) 

v e (y)<Ke- 2 m°(A(z,ej). (5.6) 

Let T a (h,z) be the cone with opening 2a, vertex at z and axis of symmetry in the h 
direction opening in the h direction. Since D is Lipschitz we can find an a and a unit vector 
h such that T a (h,z) is an exterior cone to D. Let 7 = 1/ since, and set a f = z + e^h. Then 
for e sufficiently small, one has 

B(a e ,e) C T a (h,z) n D c 
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B(z,2e) C B(a e ,(j + 2)e) 
B(a e} (7 + 3)e) DdD C A(z } (2 7 + 3)e). 

As in the proof of Theorem 3.1 of [2], it follows then that 

P y {B T2e £ dB(z,2e)) < KP y {B TD £ A(z,5e)). 

Finally (5.5) follows by the doubling lemma on harmonic measure (lemma 5.9 of [17], 
which requires d > 3). 

Next we show a similar estimate holds for a lower bound, again using the same proof as 
in [2]. By the Cauchy-Schwarz lemma we have 

v e (x) = n x (X D (A(z,e))>0) 
~ N x (X D (A(z, e)) 2 ) 

fG d (x,y)mJ>(A(z,e)ydy- 

The last line following from the Palm formula, Lemma 2.6. 
The lower estimate 

v e {x)>Ke- 2 m°(A(z,e)) (5.7) 

will follow if we show that 

G D (x, y)m°(A(z } e)fdy < Kc 2 m°(A(z } e)). (5.8) 

To do this in the d > 5 case we use two results for Lipschitz domains. First, the "3-g" 
theorem of Cranston, Fabes and Zhao (cf. [5]) shows that there exists a constant c, depending 
on the domain D } for which for all x,y, z £ D the following bound holds 

G D (x } y)G D (y } z) ( 2-d 1 I , 2 -d\ 

n <c{\x-y\ +\y-z\ ). 

Gd{x,z) 

Second, we use the following comparison between the harmonic measure and the Green 
function (cf. [17], Lemma 5.8) 

>G(y,A £ 

where M > 1 is a constant depending on D, y £ and A t is a point satisfying eM -1 < 
\A t — z\ < eM and eM -1 < d(A e} dD) (A e may be taken to be z + where T a (h,z) is an 
interior cone at z). 

First, set B = B(z } 2Mc). Then on B D D one has for e sufficiently small that 

G(x,y)/G(x,A e ) < K. 
Since m®(A(z, e)) must be < 1, it follows that 

/ G D (x, y)m°(A(z, ejfdy < KG(x, A £ ) f dy < KG(x, A e )e d < Ke 2 m°(A(z, e)). 

J BC\D JB 



M" 1 < "'''Ji < M, (5.9) 



The last inequality is a consequence of (5.9). 
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Next, 

G D (x,y)m®(A(z,ej) 2 dy 



B c nD 



<K I e d ~ 2 m D x (A(z, ej) y } G °^ A) G D (y, A t )dy (5.10) 



' B c nD X Gb{x^A ( 

<Ke d - 2 m°(A(z } e)) f (\x - y\ 2 ~ d + \y - A t \ 2 ~ d ) G D (y } A t )dy (5.11) 



B c nD 



< Ke d - 2 m°(A(z, ej) (k + [ \y - A t \ 2 ^dy 

\ JB(A e ,Mt) c 

<Ke d - 2 m°(A(z } e)) (l<+ f 



' B(A e ,Me) 

dr r 3 ~ d 

'Me 

,2„D 



< Ke 2 m^(A(z } e)). (5.12) 

Here, (5.10) holds by (5.9), and (5.11) by the "3-g" theorem. Note that (5.12) requires that 
d > 5. Therefore in that case we have by (5.7), (5.6) and the boundary Harnack principle, 
that 

v e (y) T ^ 2 m°(A(z } e)) 

— < A " rr< KK x (y,z). 

v f (x) e 2 m^(A(z,e)) 

In the case d = 4, the argument proceeds to (5.10) just as before. Choose r/ > and a 
closed circular cone C with vertex z, such that C D B(z } r/) C D c . Set D = IR 4 \ C. Then 

G D (x, y)m°(A(z, ej) 2 dy < Ke 2 m?(A(z, ej) U + f \y - A e \~ 2 G Do (y, A e ) dy) . 
B c nD v JB c r\D 7 



Thus, to show (5.8) in the case d = 4, it will suffice to show that 

\y-A £ \- 2 G Do (y,A £ )dy<K. (5.13) 



B c r\D l 



o 



We may assume that z = 0. If y,y' £ -Do and |y'| < 1, |y| > M, then Gx> (y,y') 
is dominated by a harmonic function of the form |y| _Qf 0(|y| _1 y). Separating variables in 
Laplace's equation gives that a = 1 + \A + fii f° r some (3 > depending on the opening 
angle of C. Thus a > 2, and so by the Brownian scaling and a change of variable, 

/ \y-A e \- 2 G Do (y,A e )dy = (e- 2 ) 2 |y - e" 1 A t \~ 2 G Do {y A t )e A dy 

JB c C\D J B(0,M) c nD 

<K f ^Gn^y.e-'A^dy 

J B(0,M) c nD 



oo 

2-«3 



< K / r~ 2 - a r 6 dr < oo. 
Jm 

This shows (5.13). 
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Suppose now the lemma holds for all proper subsets of A. Set a = 4|A|e. From Theo- 
rem 5.2 we have 



-Avf = kfvf -2 



\ BnC U, B v c 



2 

B,CCA 
<b,A^B,C 
BUC=A 

with kf = 2(2uf + (-l)\ A \ v f) > 0. Set 

E = {(B, C) : 0, A ^ B, C; B U C = A; \B n C\ is even}, 
O = {(B, C) : 0, A ± B, C; B U C = A; \B n C\ is odd}. 

By the Feynman-Kac formula we have 

vf(y) = E y (v?(B Ta )exp - jH kf(B r )dr] 

+ 2 J2 E y (^ Ta (Af)Wexp(- f kf{B r )dr)dt 



o 



2 J2 E y ( r( v ? v ?)( B t) e M- fkf{B r )dr)dt 
n \Jo Jo 



o 

'■Ta 



<E y (v?(B Ta )exp(- / k?(B r )d 



'o 



tr 



+ 2j2E y ( [ Ta (v B vf')(B t )exp(- f k?(B r )dr)dt 

T? \Jo JO 



E 

= 1 + 11. 

First we analyze the harmonic term /. Since kf > we have by (5.6) that 

/ <Y,E y {vt{B ra )l{B ra e dB(z t} a))) 

ieA dB(zi,a) 

In particular we have 



< 



U 3 eA-W ft A \ IWW(*0 y^?(A(^a)) 

< KY^e 2 K x {y, Zi ), (5.14) 

ieA 

since the first factors are bounded by induction. We see that the contribution from the 
harmonic term vanishes in the limit as e — > when \A\ > 2. 

To bound II we first remark that vfv^ < vf\ c v^ (or vfv^\ B if B \ C is empty), and so 
it suffices to assume that we can bound the terms of II when B D C = 0. The induction 
hypothesis applies to proper subsets of B C A as a > 4|A|e > 4|_B|e. For B D C = 0, 
B U C = A, B ^ 0, C ^ 0, we have that 
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E y [ r a (v?v?)(B t )exp(- fkf{B r )dr)dt 



< Km v[(x))E y [ r dt K?(B t , k ?( b *> z ^ ) ■ 

ieA V B C / 

As II is bounded by a finite constant times terms as above, it suffices to consider terms of 
the following type with i ^ j 

E y dtK^B^zAK^B^z^ = Gd ('■> z i)Kx ('■> z i)^j (y)- (5-15) 

Since Zi and Zj are separated by at least 8 we have that there exists a constant K depending 
only on 8 and D such that 

The "3-g" theorem gives 

n (t ,d ( yw n f r G D (y,w)G D (w,z) 

GD{K{-,Zi)){y)= hm -— dw 

Jd z ^ z < G d (x,z) 

< K f hm G ^ D{y ' Z \ (\ y - w f- d + \ w - z f- d ) dw 
J D *^*i G D (x,z) 

<KK?(y,Zi). (5.16) 

Putting together lines (5.16), (5.15), and (5.14) we get the necessary bounds for the proof 
of the lemma. □ 

Next, to finish the proof of Theorem 5.3, we again use induction on the size of A. First 
consider the case where we hit a single point. Let A = {i}. Let 8 > 0, and recall that 
t$ = rf = t d a. Where convenient, we will also write D$ for Df. Let P yz be the law of 



Brownian motion conditioned to leave a domain at z starting from y. 
We start with the following lemma 

Lemma 5.5. Let y £ Df. Then uniformly in z £ 3D 

limP^(exp(- ^ ' kf(Bt)dt)) = 1. 



Proof. We have v[{x) — > 0, so by Lemma 5.4 it suffices to prove that 
uniformly in z. But 



limP^(exp(-Ar 5 )) = l, (5.17) 

A-»0 



sup P yz (rs) < oo 

yED s ,zEdD 

since D is Lipschitz [5]. Line (5.17) then follows immediately by a well-known argument of 
Khasminsky. See Lemma 3.7 of [4]. □ 

Let x, y £ D$ 0} by the Feynman-Kac formula, for each fixed 8 < 8 we have 

v[(y) _ E y (vl(B TS )e X p-f TS 2k?(B r )dr) 

vl(x) E x (vl(B TS ) exp - 2kf(B r )dr) 
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= n aDs ^(exp - J? k?(B r )dr)K?'(y, z)vj(z)m^ (dz) 

\ J d D S v i( z ) m * s ( dz ) 

x ( f dDs <(z)™°Hdz) \ 

X \J aDi P **(eW-j; s WBr)dr)vi(z)m°'(dz)) ' 

The measure 



X e ,s(x,dz) = /, * X D G Mt(dD 8 



f aDg vi(z)m x °(dz) 

Since the boundary of Df is compact, by Prohorov's theorem any sequence e 3 has a subse- 
quence, again written ej, for which =>■ \s G M.\(dDf). Also, K® s (y,z) is continuous 
and bounded in z, for z £ Dfl dDf, when x,y G -D^. Since P^(exp — J^* kf(B r )dr) — > 1 
uniformly as e — > by Lemma 5.5, we have that for x, y G -D^ and for all 8 < 8 

lim = = / K?'(y,z)\ s (x,dz). 



This limit is then harmonic in y for y G Dk . By a diagonalization argument, we can 
assume there exists a convergent subsequence of our sequence such that the convergence 
holds simultaneously for a sequence of 8j's which converge to 0. By Lemma 5.4 we see then 
that the limit is harmonic in y with boundary value on 3D H dDf for all 8 > 0, and is 1 

at y = x. In other words, linij ^^ 3 \y) is the Martin kernel for Brownian motion in D. Thus 
all subsequences have a subsequence which converges to the Martin kernel, and so the limit 
itself exists. 

To prove the induction step, fix A and assume the result is true for all proper subsets of 
A. Therefore, if B } C C A and B, (7 / A we have 

Again by the Feynman-Kac formula and Theorem 5.2, we have a 8 such that if 8 < 8 and 
x, y G D$ then 



Consider the right-hand side as e — > 0. By Lemmas 5.4 and 5.5, the dominated convergence 
theorem applies, and thus by the induction hypothesis, 



^° Ua v l\ x ) 



BnC=9 
BUC=A 



For the upper bound, set a = 12|A|e. Then as in the proof of Lemma 5.4, for y G the 
differential equation for vf again yields the bound 

vf(y) < E y (v?(B Ta )) | 2 ^ E y (j^(vM)(Bt)dt) 



v:\x 
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As in (5.14), the first term goes to 0. Therefore we can apply the dominated convergence 
theorem to the second term as before, and get by induction that 



BnC= 
BUC=A 
B,C+% 



□ 



Theorem 5.6. Let D satisfy the conditions of Theorem 5.3, $fc £ Tk and fix x £ D. Then 
for y £ D 



1 

limN^S* | 1[(X d ,1a0 > 0) = ^— N y ($ k Mi 



k )■> 



8 = 1 

where 

aeV{N) ceo- 

Remark 5.7. By Theorem 5.3, the probabilistic representation of section 4 applies to the 
limiting measures 

x \y ) 

We take g = 0. The backbone is a tree connecting y to the n points {z i} . . . z n } } which throws 
off mass that evolves as an unconditioned superprocess. 

Proof. The idea is to use the special Markov property of the Exit measure to get a formal 
expression for the limit and then show that the convergence indeed holds. Set = 

exp ([-iy c \X k } vf}^j - 1, so that Wf = (-1) |C| \W?\. First we consider 



,8 = 1 



(l[(X D ,l A i)>0\r k j 

limN, (ff(l-exp(-A(X D ,l A8 ))) \ T k 



A— >-oo 

\i=l / 

A limN,(l+ ]T (-1)1*1 e xp-A<X D ,]Tl A ,) \T k ) (5.18) 

Q^BCN ieB 

liml+ (-l) |B| exp-(X fc ,N.(l - exp -A(X D , £ 1 A » ))) (5.19) 

°° ®^BCN ieB 

1+ £ (-l) |B| exp-(X fc ,uf) (5.20) 
1+ £ (-1)1*1 exp- £ (-l)l^(^,,f) (5.21) 

Q^BCN Q^CCB 

1+ £ (_1)|B|( 1+ 2 ^ 1 £ JJ^) (522) 

tt^BCN j=l Ci^C 2 ^-^Cj j=1 
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= 1+ E E E II^( E (-i) |sr 

tl^BCN j=l G 1 <C 2 <-<G 3 i=l CiU-uCjCBCAT 

0^C 8 CJV Vi 

2^1-1 i 

= E E (-irn^ ( 5 - 23 ) 

j=l Ci^C 2 ^!-^Cj «=1 
C 1 U-UC j =N 

21^1-1 / i \ i 

= E E (-ir Ilc- 1 ) 1 ^ 1 ni^ Ci i 

j=l c l <C 2 <-<C J \i=l / 1=1 

C 1 U-l)C j =N 
Q^C, Vi 

2^-1 / i \ 

= e e ni^i (- i ) n+EtiiCii - 

j=l Ci^C 2 ^!-^Cj \i=l / 
CiU-UCj-AT 
0^C* 8 Vi 

In the above, (5.18) and (5.22) follow from Lemma 2.2, (5.19) from Lemma 2.5, (5.20) from 
Lemma 2.3, (5.21) from Lemma 5.1, and (5.23) from Lemma 2.1. 
Thus, we need to show that the following sorts of limits exist. 

li ^^n-l'f ' (5.24, 

IL=i v i\ x ) 

To that end, for fixed k we note that for C ^ N there exists a constant K < oo for which 
v ?(')/ Yliec v l( x ) < K m Dk- Since 1 — e~ x < e x — 1 < xe x for x > we have 

JfL < - 1 < A-(A-M>exp ((.V',l)A-n»i(= 

Moreover, by Theorem 5.3 (with dominated convergence verified by Lemma 5.4) 



lim J^VfL = (X k ^ 

so also 

LL e c W) 

Set ra™ = supp,-^ riiec ^el^)' so m e — ^ 0, as e — > 0. This gives that for bounded G Fk-, 



nr=i 



< 



c/^(X fc , 1)" exp«X fc , 1)^) n^^))"'" 1 ' 



where a 8 - is the number of times i appears in the {Ck}- By Lemma 2.7, and the fact that 
x n e x < K(e 2x — 1) for x > 0, this bound is in £ 1 (N a; ) and it decreases as e — > 0. Thus we 
may apply the dominated convergence theorem to calculate limits like (5.24). The limit is 
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unless all the a 8 - = 1, so if {Ci} 3 i=1 are distinct, nonempty and U 3 i=1 Ci = N } then 



K I r 2 ~ d dr 



£ 1]" ^'0*0 " /J 1 ^ disjoint}- 

The formula for follows immediately, using Lemma 3.4. □ 

Remark 5.8. Recall that Theorem 3.1 includes the hypothesis that v N < oo. With z±, . . . , z n 
distinct and D Lipschitz, this condition follows from Lemma 5.4. It is worth observing that 
it may fail if the Z{ are not distinct. For example, if D = -6(0, 1) and z £ dD, we take n = 2, 
z\ = z 2 = z, and v 1 = v 2 = Kq(-, z). Then for C a suitable cone, with vertex at z, 

v N (x) = 2 [ G D (x,y)K°(y,z) 2 dy 

J D 

> K [ d(y, dD)[d(y, dD)'^] 2 dy 

Jcr\B(z,i/2) 

>K f" ! r i-2(d-i) r d-i dr 
Jo 

•1/2 

r 2 ~ d dr = oo, 

'o 

for d > 3. 

Remark 5.9. What if, instead of using the excursion measure N y} we use the probability 
measure P^, under which the superprocess X t has X = /i? In this case, one expects the 
form of the answer to change, as it is no longer necessary that the backbone originate with 
a single path. In genealogical terms, there may be more than one ancestor for the particles 
reaching the set {z i} . . . , z n }. 

In fact, the argument of Theorem 5.6 still applies (assuming for simplicity that ji has 
compact support in D), so that for £ jF fc , 

limP^S* | l[(X D , l A j> > 0) = ¥ -^F,($ k M? ), (5.25) 

where = ~Y^aeT(N) Wc^ai^i ^x) as before. For simplicity, drop the subscript x. It 

remains to evaluate the expressions in (5.25). 
Recall that 

P^exp -(X\ cf>)) = exp -(fx, N.(l - ej)>. 
Let a £ V(N), and take s = \a\. As in Lemma 2.6, 

8 = 1 8 = 1 

' ^ ^ exp-( / i,N.(l- e fc (</> + VA.^W))) 



o?Ai . . . o?A s " ' v v 1 ^-^ 

8 = 1 

v( e 5) e n<^N.(ej n^*.*^)))- 

7G-P({1,. ..,«}) i=l j£7(0 
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Thus 



v«)=p,(ej) e e n^^-^n^'^))) 



ieV(N) <TieV(-y(i)) i=l 

;=l,..,l7l 



e n^. N -( e J e n<**>^»> 

E f[^,N.(eJ^W)> 

-YEV(N) i=l 

V( e 5) E n<*VM; fl (ej)>, 

7£-P(iV)-Be7 



where Mj(eJ) = ^N^ejMf ). In particular 

rJV 



-ieV(N)Bei 



Thus, the conditioned exit measure X k , with Laplace functional 



1 



can be realized as a sum of two independent components. The first is a copy of the uncondi- 
tioned exit measure, which we will denote Xq. To construct the other component, we choose 
a random 7 £ V(N) with probability proportional to Ylse-y^^ X & B )- For each B £ 7, we 
then independently choose a starting point x with law — ^ey^^/i, and then independently 
generate exit measures Xg with Laplace functionals 



We then have that 



X k =X k +Y, X B- 

In other words, the backbone now consists of a branching forest, each tree of which has the 
same description as before, except that it targets a given subset of N. 



Remark 5.10. One way to generalize 3.3 would be to seek martingales of the form 



fc,(g)j 



(5.26) 



where X k, ® J = X k (E)- ■ -®X k (with j factors), and Uj(xi, . . . ,Xj) > is a symmetric function 
of j variables. Taking u n+ \ = 0, the required conditions are that 

]-^ Xt) u 3 (xi, ...,Xj) = -2(j + l)u j+1 (xi, Xi, x 2 ,..., Xj), 
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for each j. An example would be the symmetrization of the function 



•j 



aeV(N),\a\=3 i=l 



Note that the leading term u n is n-harmonic, in the sense that it is harmonic in each variable 
separately. Rather than trying to derive an analogue of Theorem 4.2 in this context, we 
simply note that there is an integral representation of the n-harmonic functions, in terms of 
products of the Martin kernel, of the form h(xi, . . . , Xj) = Yli=i K D ( x ji z j) ( see [3])- Thus a 
transform by of the above form will typically be a superposition of transforms of the type 
arising in Theorem 5.6. See [6] or [25] for other connections between n-harmonic functions 
and conditionings. 

6. Conditioning on the Brownian snake exiting at exactly n points 

If we designate n small balls on the boundary of D and this time condition the exit measure 
to charge each one, but not to charge the complement of their union, then the following formal 
calculations suggest a representation of the exit measures as before. The new process will 
again correspond to a branching process with immigration, however the immigrated mass 
is conditioned to die before reaching the boundary of D. The importance of this type of 
conditioning is that these conditionings should be minimal, in some sense, and so should 
correspond to points of some generalized Martin boundary. Recall that Af = Uig^A*. The 
analogue of Theorem 5.2 is: 



Theorem 6.1. For A C {1, 2 . . . , n}. Define 



u(x) = N x (hit 3D) = N X ((X D , 1) > 0), 
uf(x) = N x (hit Af, miss dD \ Af) 



= N X ((X D , 1 a a) > 0, (X D , l dDXA A) = 0), 
vf{x) = N x (hit each A*, miss 3D \ Af ) 



= N X ((X D , 1 A » ) > Vz G A, (X D , 1 9DW = 0). 



Then vf is twice differentiable in D and satisfies 




BUC=A 



BUC=A 



with boundary condition 




zedD\ A\ 

z e A\ 




A\ > 2, 



Conjecture . We conjecture, in analogy to Theorem 5.3, that the limit 
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exists, where 



Uv A = 4uTi-2 ]T r?rf fl , (6.2) 



BCA 



and that moreover, as in (3.2), 

= 2 ]T U^(T B rf fl ), for |A| > 1. (6.3) 



BCA 
<b,A^B 



Given this, let Mf and be as in (3.3) and (3.4), with g(x) = u(x) = N^(hit dD) and 
v % = rjc\ We will give a formal argument, along the lines of Theorem 5.6, that 



limN,($ fc | hit each A*, miss dD \ Af ) = ^ —N y ($ k MZ), (6.4) 



1 

r45"- » — iffy) 

for $/; G Tk- Note that the particle representation of Theorem 4.2 would automatically 
apply. 

For ease of notation, we will drop the e subscripts. We introduce the notation 
U A = {(X D , 1 a a) > 0, (X D , l dDXA A) = 0}, 
U A = {(X D ,l dDXA A)>0}, 

V A = {(X D , l Al > > Vz G A, (X D , l dD \ A A) = 0}, 
u A (x) = N X (U A ), 
so that u A (x) = N X (U A ) and v A (x) = N X (V A ). 

Lemma 6.2. We have the following relationships, 

(a) U A = [J^^bca V B as a disjoint union, so u A (x) = YIq^bca v B (x) . 

(b) \jja = \u — \u A , where U = {hit 3D}. 

^i^ = E ^(-i) H+|B| i^- 

Proof. To prove the third statement we rearrange terms as follows 

£ (-1)14+1*11^ = E ^ ^ (°)) 

Q^BCA Q^BCA Q^CCB 

= E ivc(-i)'^ E 

%+GCA CCBCA 

= ]r lvc(-iy A \-l)Mlc=A = l V A- 

Q^CCA 

Proof, (of Theorem 6.1) 

By (a) and (b) of Lemma 6.2 we have 

v *= ]T(_i)H+iV = ]T (-i) lAl+m (u-u B ). 

Q^BCA Q^BCA 

Thus 

Av A = Y (-l) H+|B| (Au- A« B ) 
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□ 



[ u 2 — u 2 B ) 



(-1)^+1*14(1 

Q^BCA 

4 {-l)\ A \+\ B \ {u 2 - (u - u B f) (6.5) 



Q^BCA 



4 ^ (_i)N+|B| ( 2uu 



0^£CA 

4(2u ]T (-1)1^+^1^- £ ]T ^) 2 ) (6.6) 

a(2uv a - £ (-l) lAl+lBl E v °' vD 

Q^BCA 9^C,DCB 

a(2uv a - Y E 

0^C*,DCA CUDCBC4 

4(2 ra A - E (6.7) 

0^C*,D 



4(2u^-2^ J] u C + (^)2_ ^ 



C D 
V V 



%±CdA CUD=A 
%,A+C,D 



4(2{u - u A ) + v A )v A - 4 Y 



C D 
V V . 

C'UD=A 
%,A+C,D 

Here, (6.5) and (6.6) follow by (b) and (a) of Lemma 6.2, and (6.7) follows by Lemma 2.1. 
As for the boundary conditions, we note that by inclusion 

N^hit each A*, miss dD\Af)< N^hit A^') 

for each j. If \A\ > 2, then each z G 3D the latter tends to for some j. The case A = {i} 
follows as before. □ 

Now, assuming the conjecture, we sketch an informal argument for (6.4). 

= N x ($ k E (-l) n+ Wl uB (6.8) 

Q^BCN 

= (-l) n K*(*k E (6-9) 

Q^BCN 

= ]im(-l) n N x ($ k E (-l) |B| (exp(-A(X D ,l))-exp(-A(X D ,l A? ))' 

0^BCJV 

= (-l)"N iC ($ fc E (-l) |B| (exp-(X fc ,u B )-exp-(X fc ,u))) (6.10) 

Q^BCN 



l) n n x U k exp(-(X k ,u)) E (-l) |B| (exp(X fc ,u B )-l) 



0^BCJV 
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l) n N x (^ k exp(-(X k ,u))(l+ ]T II ex P( X 



= (-iyn x (<$> k eM-{x\u))(i+ ]T II (i + (ex P (x fc ,^)-i; 

Set IT^ = exp(X fc , u c ) - 1, so by Lemma 2.2 

n (i+w g )=i+y, e ri^- 

tl^CCB j>l C\-<--<Cj i=l 

0^C* 8 C£ Vi 

Again, we note that ^0^bcjv( — = — 1 by Lemma 2.1. Thus we may continue the above 
calculation to arrive at 

0^£CiV j>l Ci^-^Cj «=1 

0^C* 8 C£ Vi 



irN,($ fc exp(-(X fc ,u))(E E E ("^ 



\B\ 



j>l Ci<-<C 3 8=1 UC.CBCN 



j>l Cl^-^Cj 8 = 1 

0^C* 8 Vi 
UC 8 =JV 

the last equation being a consequence of Lemma 2.1. Therefore, if we assume there are no 
problems with the convergence, we have that 



e^O e^O V N (x) 



j>l C l <---<C ] i=l LLN eX ' N 

UC,=N 

N,($* exp(-(X fc , «» Z,eV(N) IW<**, if » 



r^(x). 

The last line follows as the terms with the C's not disjoint disappear in the limit. 



7. Historical processes 

In previous sections we have, for simplicity, avoided the use of historical processes, even 
though this would have given us better results. In this section we remedy this problem, and 
sketch how those earlier results can be strengthened. 

We start by adding to the material of section 2.4. Again, let D C IR d be a domain, and 
take bounded, smooth domains D k ff D, with X k the exit measure from D kl and Tfc(£) the 
exit time of £ from D k . Recall that (W S} ( S ) denotes the Brownian snake. Let W\ denote 
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the stopped snake, W*(it) = W s (t A it). Define the historical process, of paths killed upon 
exiting D k} to be 



(#*,$> = j $(W:)l Tk < t dL t (s) 



where now $ is a measurable function on path space. 

We have repeatedly used the cr-fields Th of information determined by the super Brownian 
paths prior to exiting D k . Their formal definition is that 

T k = a{H k | t > 0}. 

Let 

poo 

^($.) = exp- / (H k ,$ t )dt, (7.1) 



where t \— > $ t is positive and continuous. Since H k is continuous under N^, or under any law 
absolutely continuous with respect to N X} such laws are characterized on T k by the expec- 
tations of random variables /jq($.). In particular, our approach to extending Theorem 4.2 
will be to take the calculation of the e k ((f)) = exp — (X k ,(f)) under M^, and extend it to the 

To do so, we need to record versions of the Palm formulae in this context. A first step 
towards doing so is to redefine our underlying measures, allowing them to start from a path 
rather than a single point. If w is a path, and t > 0, we let ~N Wjt be the excursion measure of 
the Brownian snake, started from the stopped path w t (s) = w(s A t). In other words, under 
N m t, we have that Wq = w f and that the law of (. is that of a Brownian excursion above 



w,t 



level t. Thus L u = for u < t, so that H k = for u < t. Alternatively, we can obtain N. 
directly as the image of ~N W ^ under the mapping £ i— y £, where 

yt(s-t), s>t. 

We write P w>t for the probability under which B s = w(s) for s < t, and for which B t +. has 
the law of a Brownian motion started from wit). For s < t, let 

ft 

Jz( ±\ l,k{ 



Af Stt (e ((f)), h k ($.)) = exp - j du 4N B , U (1 - e k (<f>)h k u ($.)). 

The basic Palm formula (2.3) then takes the form 

n w , t ((X k ,(f))e k (4,)h k (M>.)) = E w , t ^(B Tk )M t , Tk (e k (^),h k (^.))), (7.2) 

whenever T k (w) > t. This may be proved as in Proposition 4.1 of [20], using the point of 
view of chapter 4 of [8]. 

Using (7.2) in place of (2.3), the following extended Palm formula may be proved, exactly 
as in Lemma 2.6. 

Lemma 7.1. Let N = {1,2,... n} , n > 2. Let {tpi} be a family of measurable functions, 
and suppose that T k (w) > s. Then 

ieN 
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l - ]T E WtS (4 r dtM Stt (e k (<f>),h k ($.)) 

71 /r r~ at ^ S 



MCN 
0,JV^M 

ieM ieN\M 

We now define a probability Q w>s under which T t is a measure-valued process living on 
a branching tree. The tree starts as a single particle following w, but evolves as before 
following time s. In other words, T t = S w ^ for t < s, and T s+ . has the same law under Q w>s 
as T. under Q w ( s ). Of course, we actually need a historical version of this construction, so 
we let T t be the atomic measure on path space, which gives unit mass to each path which 
is the history of an atom of T t . In the notation of (4.2), 

nt 

T t (du) = y^S x i(du). 

i=l 

Tki 



To define N W)S , we again form a Poisson random measure N (o?x,c//i) with intensity 

dt I 4T k (du)N UJ 4H k G d X ,X k G dp), 



oo 



where x belongs to the space of paths of measure- valued processes, and ji belongs to the space 
of measures on dDk- We then realize the exit measure under N WyS as X k = j ^N k (dx,d^), 
and the historical superprocess as H k = j xN k (dx,d/j). In concrete terms, 

n w ,s(e k (<P)h k s (^)) = Q w ,s(exp- j ^4(f*,N,t(l-e fc (^($)))). (7.3) 

•J S 

The generalization of Theorem 4.2 is then 
Theorem 7.2. In the notation of Theorem 4-2, if v N < oo then = on Tk- 
Proof. The result to be shown is that if Tk(w) > s, then 

M^e^)/^)) = N WiS (e k (<i>)h k ($)), 

where 

ML, s ( e fc ^)^($)) = 1 N w , s ( e fc (^)^($)Mf). 

[W{S)) 

In fact, the argument proceeds exactly as before, with induction being used to prove a 

generalization of (4.6). The necessity of considering and ~N W>S , rather than just and 

Re, arises from the induction. □ 

Similar arguments can be used to show that other measure equations, proved before simply 
for expectations of random variables e^, can be extended to arbitrary elements of Tk- We 
leave the details to the interested reader. 
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